Differentiating the incidence and burden of HIV by age among women who sell sex: a systematic review and meta‐analysis

Abstract Introduction Young women who sell sex (YWSS) are at heightened risk of HIV acquisition and transmission and are among the least engaged in HIV services. There is insufficient epidemiologic evidence characterizing the burden of HIV among YWSS, particularly as compared to older WSS. These data are needed to design and tailor effective HIV prevention and treatment programmes for this population. Methods We conducted two parallel systematic reviews and meta‐analyses to define both the immediate and long‐term HIV risks for YWSS, including among women engaged in sex work, survival sex and transactional sex. In the first review, we identified and synthesized published studies of HIV incidence comparing estimates for cisgender women ≤24 years of age versus >24. In the second review, we identified and synthesized studies of HIV prevalence, comparing estimates for cisgender women who initiated selling sex <18 versus ≥18 years. In both reviews, we completed a search of four databases for articles in any language and any geographic area published from 1 January 1980 until 12 February 2021. Included articles were assessed for quality and a random effects model was used to calculate pooled effect estimates for each review. Results and Discussion We identified 12 studies for the HIV incidence review and 18 studies for the HIV prevalence review. In a meta‐analysis, HIV incidence was elevated in younger (5.3 per 100 person‐years [PY]; 95% confidence interval [CI]: 3.5, 7.1) compared to older women (2.8 per 100 PY; 95% CI: 1.7, 3.9), although CIs overlapped. HIV prevalence among those who initiated selling sex <18 years of age (28.8; 95% CI: 18.9, 38.7) was higher than those who initiated later (20.5; 95% CI: 12.4, 28.6). Conclusions These companion reviews offer an important perspective on the relative HIV risk of engaging in selling sex at a younger age. Our findings highlight the unique and intersectional challenges YWSS face, and the importance of ensuring that health services are tailored to meet their specific needs. Research and programming should routinely stratify data into meaningful age bands to differentiate and intervene within this population.


I N T R O D U C T I O N
Due to the increased availability of biomedical technologies for the prevention and treatment of HIV, there have been significant decreases in HIV incidence, morbidity and mortality observed globally [1,2]. Despite these successes, epidemiologic gains have not been evenly distributed across populations. As a result, the global HIV agenda has seen an increased focus on improving HIV prevention and treatment for key and priority populations who remain at heightened risk of HIV acquisition and transmission and are marginal-ized from existing services, including both young people and women who sell sex (WSS). Key populations, including sex workers, men who have sex with men, transgender persons and people who use drugs, have among the highest risk of contracting and transmitting HIV irrespective of the epidemic type or local context and frequently lack adequate access to services [3]. Other vulnerable populations who have also been prioritized by the global HIV response include adolescent girls and young women ages 15-24, who are particularly vulnerable to HIV infection in Southern and Eastern Africa [3]. While young women and WSS are often distinguished as distinct populations in both programmatic and research settings, there are substantial overlapping vulnerabilities. Between 20% and 40% of WSS first initiate sex work during adolescence and a review of HIV incidence in young women in 10 countries found that some of the highest incidence rates were among female sex workers aged 18-24 years in South Africa (13.20 per 100 person-years [PY]) and in Zimbabwe (10.80 per 100 PY) [4,5]. Therefore, addressing the intersectional needs of young women who sell sex (YWSS) is critical for optimizing HIV prevention and treatment, meeting UNAIDS 2025 targets to ensure that 95% of women and girls of reproductive age have their HIV and sexual and reproductive healthcare service needs met, and advancing the global agenda of ending HIV/AIDS [6,7]. Achieving these goals will also require efforts that address the underlying social and structural determinants that drive HIV acquisition and transmission for this population [8,9].
Evidence quantifying the burden of HIV among YWSS across settings is limited, and few published studies explore the longer-term epidemiologic implications for women who initiate selling sex at a young age, precluding a full understanding of HIV incidence and prevalence for WSS during adolescence and early adulthood. Moreover, differences in HIV risk may vary by types of sexual exchange, including formal sex work, survival sex or transactional sex, each of which typically reflects different motivations, partner power dynamics and economic circumstances [10][11][12]. Finally, young women under the age of 18 cannot be legally defined as sex workers and, therefore, are largely excluded from epidemiologic research due to ethical concerns and legal constraints [4,13]. These gaps in evidence challenge effective resource allocation and service delivery, and data are urgently needed to better assess and meet the otherwise unmet HIV prevention and treatment needs of YWSS.
The objective of this research was to define both the immediate and long-term HIV-related risks for YWSS. To do so, we conducted two systematic reviews in parallel, responding to two related research questions. The first review and meta-analysis aimed to assess whether YWSS have a higher HIV incidence as compared to older WSS. The second review aimed to determine if young age at initiation of sex work was associated with a higher HIV prevalence as compared to those who initiated sex work at an older age. The findings of the two reviews are presented and synthesized together to provide complementary insights into HIV risks related to young age among WSS.

M E T H O D S
We conducted two independent systematic reviews utilizing the standard protocol for Preferred Reporting Items for Systematic Reviews and Meta-Analysis (PRISMA; see Figures 1 and 2 for the screening flow diagram and Appendix S2 for the PRISMA checklist) [14]. Reviews were conducted concurrently. In the first review, we identified and synthesized all published studies of HIV incidence among WSS that reported disaggregated estimates for women ≤24 years of age (hereafter referred to as the "HIV incidence review"). The categorization of young women under 25 years was selected to align with the UNAIDS definition of adolescent girls and young women which is 15-24 years of age [1]. In the second review, we identified and synthesized all published studies of HIV prevalence among WSS that reported disaggregated estimates for women who initiated selling sex ≤18 years of age (hereafter referred to as the "HIV prevalence review"). We used a similar methodological approach for each review and present our modifications to eligibility criteria and search strategy which differentiate each research question. Metaanalytic methods were used for both reviews.

Eligibility criteria
Studies published from all countries were eligible for inclusion without restriction by geographic region. Eligible studies were those that (1) included cisgender women who sell sex, including sex work (street-based, venue-based, those working at sex work establishments), cisgender women engaged in survival sex work, and cisgender women and girls under the age of 18 engaged in transactional sex; (2) reported sufficient detail on sampling methodology, study design, methods of evaluation and analysis of results; (3) included ≥20 individuals who identified as women or girls and were not limited to men or transgender populations; and (4) were published between 1 January 1980 and 12 February 2021. Systematic reviews and conference abstracts were excluded. While survival sex, sex work and transactional sex represent different types of commodified sex, they have all been associated with increased risk of HIV and are motivated by different factors which can influence risk [10,11,[15][16][17][18]. Further, young women ages 18 years and older are considered adults and, therefore, fully able to determine their engagement in the sale of sex; however, the sale of sex by any person under the age of 18 is defined under international law as sexual exploitation [13]. Definitions of YWSS (e.g. sex work, women engaged in survival sex work and transactional sex) were based on the definitions used by the authors of the studies and were all included to represent a fairly wide definition of any type of sexual exchange that could include YWSS over time. Transactional sex was defined based on studies that used this term (typically as the exchange of sex for money, drugs or other gifts) and was included because transactional sex has been associated with a higher risk of HIV among young women and is often used as a proxy for selling sex for those legally under age 18 [11,13,19]. Survival sex was defined based on the use of this term in the study of interest.

Methodology unique to HIV incidence review
In addition to meeting the aforementioned eligibility criteria, studies were required to report biologically measured estimates of incident HIV infection. These included both prospective and retrospective longitudinal cohort studies as well as cross-sectional studies where HIV incidence was captured using an incidence assay (e.g. BED capture enzyme immunoassay). Randomized trials reporting the effects of intervention were included if data were disaggregated by the study arm; in this case, we reported incidence data from the control arm only. Additionally, studies were required to include young women ≤24 years of age and present results stratified by age. We extracted available quantitative data related to HIV incidence for each study. Key variables were extracted overall and by age strata, and included sample size, person time contributed, number of HIV seroconversions and HIV incidence per 100 PY. For studies where person time was not reported, we extracted risk ratios as reported in the study.

Methodology unique to HIV prevalence review
Studies eligible for inclusion were those that met the eligibility criteria above and reported HIV prevalence, as measured by self-report or through biological testing. All quantitative study designs were eligible. For randomized trials testing interventions, we extracted pre-intervention measures of HIV prevalence. Studies were also required to report the age at which selling sex was initiated, include women who initiated selling sex at ≤18 years of age and to present data stratified by age at initiation of selling sex. We extracted all available quantitative data related to HIV prevalence, including sample size, the number living with HIV, HIV prevalence and HIV prevalence odds ratio.

Study selection
A search of Pubmed, EMBASE, Web of Science and Sociological Abstracts databases was conducted for articles in any language and in any geographic area published from 1 January 1980 until 12 February 2021. Abstracts not in English were translated via Google Translate [20]. These databases were selected a priori as those most likely to contain relevant citations. The search strategies were piloted to assess whether pre-identified relevant studies had been captured and were refined accordingly. The final search strategies are described in detail in Appendix S3. Abstracts were imported from each database, duplicates removed and then screened for eligibility based on title and abstract by two independent reviewers using Covidence online software [21]. Articles flagged for inclusion by one or both reviewers underwent a full-text review. The final decision for inclusion was based on an agreement between the two reviewers, with a third reviewer available in the event of a discrepant decision by the first two reviewers. We also conducted a hand search of references of the citations which met inclusion criteria, to identify any additional relevant studies. We did not identify any additional articles this way, but this provided confirmation that we had included all relevant articles. Descriptive data about the study and sample were extracted, including study design, methods, type of analysis, population/inclusion criteria, the definition of selling sex, age groups reported, time period, geographic location and any other relevant outcomes. We also extracted information relevant to either HIV incidence or HIV prevalence, depending on the outcome.

Quality assessment and risk of bias
We assessed the strength of each included study using a modified quality assessment tool, the Global HIV Quality Assessment Tool for Data Generated through Non-Probability Sampling (GHQAT) [22]. The tool was specifically developed to evaluate the quality of HIV epidemiologic evidence generated using non-probability methods, which is common in research studies with key populations. Members of the Global HIV Research Group that supported the development of the tool include collaborators from the Johns Hopkins School of Public Health, The United States Agency for International Development (USAID), The Global Fund to Fight AIDS, Tuberculosis, and Malaria, and UNAIDS. Based on the National Heart, Lung, and Blood Institute's Quality Assessment Tool for Observational Cohort and Cross-Sectional Studies, the tool evaluates study design, implementation and the strength of core epidemiologic indicators common in HIV-related key populations research. Each study was scored using the GHQAT by two independent reviewers. Studies received a total score which comprised an evaluation of study design, study implementation, and indicator-specific criteria. For the HIV incidence review, measures of reported HIV incidence were also graded; for the HIV prevalence review, measures of HIV prevalence were graded. Overall scores were then assessed as good, fair, or poor based on the tool's established criteria. Studies with an overall score of "poor" were excluded from the metaanalyses in sensitivity analysis to assess the undue influence of these observations on results and findings.

Data synthesis and statistical analysis
The primary measure for the incidence review was HIV incidence per 100 PY. The primary measure for the prevalence review was HIV prevalence. Data were organized both in tabular form and graphically using forest plots. Where 95% confidence intervals (CI) were not reported, asymptotic 95% CIs were constructed for each point estimate using Rothman's EpiSheet [23]. Age-stratified estimates for each review were visually inspected to assess potential trends by age. Subanalyses were conducted to examine how estimates varied across geographical region and definition of YWSS (e.g. sex work, women engaged in survival sex work and transactional sex). For the incidence review, we conducted a meta-analysis to calculate the overall HIV incidence rate for both younger and older age groups separately. Given the heterogeneity in age groups that were reported across studies, we categorized incidence rates as occurring among YWSS less than 25 years and among those greater than or equal to 25 years. Incidence rates for multiple age groups above age 25 were combined within studies (e.g. if both 30-39 years and 40-49 years were reported in a single study), and 95% CIs were recalculated for these combined groups. Five estimates were reported for YWSS ages 20-29 years. Given that these estimates did not explicitly fit within either age category (<25 vs. ≥25), they were reported as a separate group. Two additional studies reported an estimate for ages 12-25 years and were not included in the meta-analysis because it overlapped with the cut-off (≥25). We conducted a sensitivity analysis including these two studies. We found that when recalculating the incidence estimate, the results were similar to when these studies were not included (∼0.50 increase in estimate). We used a random effects model to calculate a pooled incidence rate. A weighted average of the pooled estimates was calculated using a random intercept model to account for between-study variation [24,25].
For the prevalence review, we similarly conducted a metaanalysis, dichotomizing all point estimates of HIV prevalence as those that were captured among women who initiated selling sex at less than 18 versus 18 years or older. This cut point was chosen based on the legal age of consent and to maximize fidelity to included studies, most of which similarly dichotomized prevalence estimates at 18 years. We used a random effects model to calculate a pooled prevalence estimate, using the same approach as the incidence review [24,25]. Age was routinely reported as a categorical variable across the studies that were included in these analyses, with HIV incidence and prevalence reported for each age category. Therefore, we report age categories and not continuous measures of age throughout.

Role of the funding source
This study was funded by the National Institute of Mental Health. The study sponsor had no role in the study design, implementation or the decision to submit the paper for publication.

HIV incidence review
After the exclusion of duplicate citations, a total of 3725 studies were identified from four databases. The title and abstracts were screened for eligibility, and 239 were retained for full-text review. Reasons for exclusion are reported in Figure 1. A total of 12 studies were eligible for inclusion [26][27][28][29][30][31][32][33][34][35][36]. Only 11 studies were included in the meta-analysis because one study included age stratification which overlapped with our age cut-off. This study was not included in our primary analysis but was included in the sensitivity analysis. These studies are summarized in Table 1, grouped by geographical region. Studies were conducted in sub-Saharan Africa (N = 6) and in Asia (N = 6). Most studies were prospective (N = 9) or retrospective cohort studies (N = 1), and two were cross-sectional studies that used an incidence assay. Of the 12 studies, seven defined WSS as self-reported commercial sex work, but one of these also measured "indirect sex work" or transactional sex. Three studies used definitions of transactional sex. The last two studies were not explicit in how they defined selling sex or sexual exchange. Overall HIV incidence ranged from 1.05 per 100 PY in Thailand [28] to 10.8 per 100 PY in Zimbabwe [29]. All studies were published between 1994 and 2018, with seven studies published since 2010. Estimates from studies that used definitions of transactional sex ranged from 2.5 [27] to 6.2 [36] per 100 PY for younger groups, while estimates from studies that used definitions of commercial sex ranged from 1.0 [30] to 11.1 per 100 PY for younger age groups [33]. In studies with an older comparison group (N = 11), the incidence was higher in younger compared to older women in six studies, comparable in three studies and higher in older women compared to younger women in two studies (Table 1). Among studies in African countries with an older comparison group (N = 5), the incidence was higher in younger compared to older women in two studies, comparable in one study and higher in older women compared to younger in two studies. Among studies in Asian countries (N = 6), the incidence was higher in younger compared to older women in four studies and comparable in two studies. In a meta-analysis of pooled estimates by age group, the incidence rate for young women below the age of 25 years was 5.3 per 100 PY (95% CI: 3.5, 7.1; Table 3 and Appendix S1 Figure 1). There was substantial heterogeneity between studies (X 2 = 38.62, degrees of freedom = 8, p<0.001) with 79.3% of the variation across studies estimated to be due to heterogeneity rather than chance. In a sensitivity analysis, where we included YWSS ≤25 (two additional studies), the incidence rate for young women less than or equal to 25 years was 5.4 per 100 PY (95% CI: 3.3, 7.4; Appendix S1 Figure 2). Among women greater than or equal to 25 years of age, the incidence rate was 2.8 (95% CI: 1.7, 3.9; Appendix S1 Figure S3) and was similar in a sensitivity analysis excluding two estimates (2.8 per 100 PY; 95% CI: 1.7, 3.9; Appendix S1 Figure S4). Again, there was substantial heterogeneity between studies (X 2 = 47.02, degrees of freedom = 7, p<0.001, I 2 = 85.1%). We estimated the incidence rate for young women aged 20-29 years to be 4.4 (95% CI: 1.3, 7.4; Appendix S1 Figure S5), again with substantial heterogeneity (X 2 = 222.91, degrees of freedom = 4, p<0.001, I 2 = 98.2%).

HIV prevalence review
After the exclusion of duplicate citations, a total of 1113 studies were identified from the four databases. The title and abstracts were screened for relevance, and 70 were retained for full-text review ( Figure 2). Reasons for exclusion are reported in Figure 2. A total of 18 studies were eligible for inclusion. These studies are summarized in Table 2, grouped by geographical region [17,[37][38][39][40][41][42][43][44][45][46][47][48][49][50][51][52][53]. Most studies were in Asia (India N = 4, Nepal and Thailand N = 2), North America (Canada N = 2, the United States = 2 and Mexico), sub-Saharan Africa (Lesotho and Rwanda) or Latin America (Ecuador and Guatemala), with one study in Ukraine and one study in Iran (Table 2). Almost all studies were cross-sectional (N = 16) but one study was a randomized trial and another was a cohort study. Definitions of WSS varied; nine studies used definitions of formal sex work, one used survival sex and eight studies used definitions of transactional sex. Overall HIV prevalence ranged from 5.2% in Iran [44] to 42.9% in Rwanda [48]. Studies were published between 1995 and 2020 with 13 studies published since 2010. The prevalence of HIV was generally higher in women who initiated selling sex at an earlier age compared to those who initiated later (N = 12), although four studies found comparable percentages in older initiators in Mexico, India, Ukraine and Kenya. Only four studies reported data on the duration of selling sex by age of initiation, an important source of potential confounding. Five of eight studies using measures of transactional sex found a higher HIV prevalence with initiation at younger ages. A total of 6/9 studies with measures of self-reported sex work and 1/1 of studies with measures of survival sex reported a higher prevalence in those initiating at an earlier compared to a later age.

Quality assessment results and risk of bias
All 12 studies identified for inclusion in the HIV incidence review were scored by two independent reviewers using the GHQAT. The overall combined median score was 12 [IQR 11,13) out of a possible 16 points, and all studies were assessed as being of similarly "fair quality." For the HIV prevalence review, the overall combined median score was 10 [IQR 9, 10] out of 13 total points. With the exception of one study that was rated as "poor," all studies were judged to be of "fair" quality [40]. Results from the quality assessment are included in Appendix S4. Given the limited number of studies originally identified for inclusion, and the consistency in quality assessment scores using the GHQAT, no further changes were made to the meta-analyses based on these findings.

D I S C U S S I O N
This study is among the first to synthesize available literature on HIV incidence and prevalence in women who sell sex by age. These companion reviews identified 12 studies with estimates of HIV incidence stratified by age, and 18 studies with estimates of HIV prevalence stratified by age of initiation of selling sex. Most studies demonstrated a higher incidence of HIV among YWSS compared to older WSS, and a higher prevalence of HIV among those who began selling sex as a minor. In a meta-analysis of pooled estimates by age group, HIV incidence was elevated in younger compared to older women, although CIs overlapped. Our meta-analysis of HIV prevalence by age at initiation yielded a similar result, whereby the prevalence of HIV among those who initiated under 18 years of age was higher than those who initiated later, although CIs overlapped. Together, findings suggest that young women who sell sex may be at a higher risk of HIV than older women who sell sex and are seemingly more likely to also have a higher burden of HIV later in life.
Targeted key population programmes have historically focused on providing services for adults ages 18 years and older without a specific focus on younger adults or inclusion of those under 18, although in some settings have expanded to offer tailored HIV programming for adolescent girls and young women. YWSS face intersectional challenges and needs related to both age and being a woman who sells sex, including not only HIV prevention and treatment, but also gender-based violence, economic challenges, mental health, housing, education and parenting needs [5,29,[54][55][56]. Needs vary tremendously during the transitional period of adolescence and young adulthood, and services designed for adults will not meet all of these varied and dynamic needs. Yet, the evidence is lacking on HIV risk specifically among YWSS, and as a result, there remains a lack of clarity on how best to meet the myriad needs of YWSS, and few HIV programmes have adequately addressed this [19,57,58]. The elevated HIV incidence seen in this review/meta-analysis, coupled with a high prevalence among those who began selling sex at    an earlier age in this review/meta-analysis, reinforces that YWSS are exposed to HIV early on and there is thus a brief window of time in which to reach them with prevention services [19]. The high incidence of HIV among YWSS in this review/meta-analysis along with lower levels of access to care and viral suppression in other studies also suggests opportunities for onward transmission to clients and other partners, as the risk of onward HIV transmission is greatest for newly acquired infections [5,29,59,60]. Younger female sex workers are less likely to know they were HIV positive, to report any Antiretroviral therapy (ART) exposure and to be virally suppressed compared to older female sex workers [54,61]. Thus, minimizing the duration of time of undiagnosed and untreated HIV infection among YWSS is critical to meeting the global HIV agenda of HIV elimination. In general, a lack of high-quality epidemiologic data for key populations challenges programmes from meeting the essential needs of sex workers, including YWSS [62,63]. Further, HIV programmes often do not disaggregate routinely collected data by age making it difficult to determine how many YWSS are accessing services [13]. Increased efforts are needed to routinely disaggregate and report data by meaningful age bands to ensure that programmes are sufficiently reached by YWSS. Overall, the quality of evidence in our review was moderate (graded as "good" or "fair" using the GHQAT) given that there were very few studies, most had smaller sample sizes within each age category, and there was substantial heterogeneity across studies.
We documented substantial heterogeneity across studies, although the small number of studies and inconsistencies in the reporting of key variables precluded a more nuanced anal-ysis of these contextual differences. Heterogeneity may be due to inherent differences between study settings, measurement of selling sex, measurement of HIV outcomes or other differences. We were unable to do any sensitivity analyses by country given the heterogeneity across studies and the small number of studies. The high HIV prevalence observed in most studies may be due to higher incidence among YWSS or a longer duration of time selling sex among those who initiated earlier. With few studies reporting duration of time selling sex, we were unable to differentiate between the two. Further, those studies that focused on formal sex work generally reported a higher prevalence of HIV overall compared to those focused on transactional sex alone. This is consistent with prior literature indicating that sexual exchange occurs along a spectrum and the motivations underlying formal sex work versus transactional sex may vary leading to differences in HIV risk [11,15,31]. There is also evidence that those who engage in transactional sex but do not self-identify as sex workers may be less likely to access HIV and other health services [64], and they may, therefore, be a population in particular need of targeted services.
There were other limitations to this systematic review and meta-analysis. First, the review identified a relatively small number of studies with data on HIV incidence disaggregated by age, and on HIV prevalence by age of initiation of selling sex. In the resulting meta-analysis, YWSS had a high pooled incidence of HIV, but CIs overlapped with those reported for older WSS. Further, due to the small sample size of included studies in the meta-analysis, we had limited precision and limited power to fully identify differences between these groups. Second, only two regions of the world were represented in the incidence review, further limiting generalizability or transportability beyond those regions. This speaks to the need for more routine collection of data to ensure we are fully able to meet the needs of this vulnerable population on a global scale. Third, the meta-analysis stratified estimates based on age cut-offs that allowed for the largest number of studies to be included as well as the standard definition used to define adolescent girls as young women. However, studies often used inconsistent age bands making it difficult to pool and compare estimates. Given the range of geographic settings, age categories and definitions of selling sex, we also did not assess meta-biases or confidence in cumulative evidence. Fourth, HIV prevalence and incidence will be influenced by HIV prevention interventions, social determinants of health and factors related to access to care that will vary by region and the year in which the study was undertaken. More data are needed on social determinants of health and access to care across studies to better understand these contextual factors. Fifth, we did not do a comprehensive review of the grey literature and may have missed information that was not published in peerreviewed journals.

C O N C L U S I O N S
These companion reviews offer an important perspective on the relative HIV risk of engaging in selling sex at a young age, as compared to older WSS. We identified an increased burden of HIV among YWSS, with generally higher HIV inci-dence compared to older WSS, and higher HIV prevalence in women who initiated selling sex at an earlier age, although CIs did overlap. Our findings highlight the unique challenges YWSS face, and the importance of ensuring that service delivery addresses their specific needs. Research and programming should routinely stratify data into meaningful age bands to differentiate and intervene within this population. Ensuring we meet their needs is a public health imperative and critical to meeting our global HIV goals.

C O M P E T I N G I N T E R E S T S
The authors declare no competing interests.

A U T H O R S ' C O N T R I B U T I O N S
MCDS, KBR and SN conceptualized and designed the study. MCDS and KBR searched and reviewed all studies for both literature reviews. CL did the metaanalysis for incidence and prevalence. MCDS led writing but all authors contributed to the writing and reviewed and edited the manuscript.

A C K N O W L E D G E M E N T S
We would like to thank Kalai Willis, Danwei Yao, Mary Anne Roach and Elizabeth Dura for their support with the quality assessment.

F U N D I N G
This study was funded by the NIH through K01MH110316, K01MH129226 and R21HD106583. CL was supported by the National Institute of Mental Health under Award Number F31MH128079.

D ATA AVA I L A B I L I T Y S TAT E M E N T
Data are contained within the manuscript or available by contacting the corresponding author.